Tuning of High-Q Photonic Crystal Microcavities by a Subwavelength Near-field 
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We demonstrate that the resonance frequencies of high-Q microcavities in two-dimensional pho- 
tonic crystal membranes can be tuned over a wide range by introducing a subwavelength dielectric 
tip into the cavity mode. Three-dimensional finite-difference time-domain (FDTD) simulations 
show that by varying the lateral and vertical positions of the tip, it is possible to tune the resonator 
frequency without lowering the quality factor. Excellent agreement with a perturbative theory is 
obtained, showing that the tuning range is limited by the ratio of the cavity mode volume to the 
effective polarizability of the nano-perturber. 

PACS numbers: 42.70.Qs, 42.50.Pq, 42.60.Da, 42.64.Pq, 42. 25. Fx. 



Solid-state optical microresonators are of great inter- 
est to a wide range of fields such as bio-sensing, non- 
linear optics, low-threshold lasers and cavity quantum 
electro-dynamics 0, Q . Two cavity properties that are 
commonly desired, but often incompatible, are very high 
quality factors (Q) and very small mode volumes. Pho- 
tonic crystal microcavities in thin semiconductor mem- 
branes ISl^iSij promise a good compromise between these 
features. One of the difficulties associated with such 
monolithic ultrasmall resonators is to match their res- 
onance frequencies with those of interest in a given ap- 
plication. This is especially a concern since fabrication 
tolerances make it nearly impossible to realize the exact 
design parameters. Temperature tuning can sometimes 
be used to meet a resonance condition. For example, 
Yoshie et al. used this technique to control the coupling 
of a quantum dot with a microcavity . This was pos- 
sible due to the temperature sensitivity of the quantum 
dot. In general, however, it is desirable to independently 
tune the cavity resonance without manipulating the sys- 
tem that it couples to. In this Letter we show that the 
introduction of a subwavelength dielectric object, such as 
a scanning probe, can achieve this. 

A Scanning Near-field Optical Microscope (SNOM) tip 
has recently been used to image the intensity distribu- 
tion in a low-Q photonic crystal microresonator 0. The 
central idea in our current work is to investigate how the 
presence of an external dielectric nanometer-sized object, 
such as a SNOM tip, can modify the resonance condition 
of a photonic crystal microresonator (see Fig.^a)). This 
scheme is analogous to a technique commonly used in 
microwave engineering whereby the frequency of a res- 
onator is adjusted by insertion of dielectric stubs 
In that case, because microwave resonators are typically 
closed, scattering from the object does not cause any loss, 
and the cavity Q is determined by absorption. However, 
when a dielectric object is placed into the field of an 
open cavity, such as a photonic crystal slab resonator. 



scattering dominates. Indeed, it has been shown exper- 
imentally that coupling of a glass fiber tip can shift the 
resonance frequency of a high-Q microresonator at the 
cost of reducing its quality factor j^. Since a photonic 
crystal microcavity has a very small mode volume and 
relies on the precise arrangement of dielectric material at 
the subwavelength scale, one might expect that introduc- 
ing the slightest external object would spoil the quality 
factor. In what follows we show that this is not neces- 
sarily the case. It turns out that a small dielectric object 
acts in the same fashion as an atom which, within clas- 
sical linear dispersion theory, changes the frequency of a 
high-finesse cavity according to its polarizability and its 
position relative to the nodes of the cavity mode 0, ^3 ■ 

We consider a dielectric cylindrical tip near a 
membrane-type photonic crystal, consisting of a high in- 
dex slab that is perforated with a hexagonal lattice (lat- 
tice constant a = 420 nm) and surrounded by air. We 
assume a slab dielectric constant e = 11.76, a thickness 
of 250 nm, and a hole radius r = 0.3a. As Fig. Qla) 
shows, the resonator is formed about a defect of reduced 
radius r = 0.15a, similar to the optimization in Ref. 
By reducing the radius of just two holes to r = 0.23a on 
either side of the defect and then shifting them outwards 
by 0.11a, we create a nondegenerate dipole mode with 
a mode profile as shown in Fig. ^b). The mode has a 
frequency Loa/2nc — 0.284 in the center of the 2D band 
gap. The Q is around 13000, corresponding to a reso- 
nance linewidth of 15 GHz at A = 1500 nm. Crystals 
over 5x5 /Ltm in lateral size surrounded by up to 1 /im 
of air were simulated using the 3D FDTD method with 
Liao's absorbing boundary conditions Compu- 
tational meshes were 14 or 20 grid points per a parallel 
to the membrane and had doubled resolution normal to 
the membrane Grid-cell volume averaging of the 

dielectric constant was employed to reduce staircasing 
errors Quality factors and cavity mode frequencies 
were obtained by fitting a damped harmonic wave to time 
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FIG. 1: (Color) (a) Schematic of the system under study. A cylindrical silicon tip with a diameter of 125 nm is placed at 
height ztip and lateral position r|| in front of a photonic crystal slab. Reduction of diameters for 3 pores (shaded red) and 
shift of two of them yields a cavity with a resonance at A = 1500 nm and a quality factor of 13000. (b) Intensity distribution 
|Eo|'^ of the unperturbed mode in a plane 30 nm above the slab. The black circles display the positions of pores (red circles 
for tuned pores) whereas the small white dots mark the original locations of the pores in the hexagonal lattice, (c) Calculated 
resonance frequency shift induced by the tip placed 30 nm above the slab, as a function of its lateral position, (d) The inverse 
Q corresponding to (c). The unperturbed cavity has 1/Q — 7.5 ■ 10~^. (e) Contribution of the ^-component |i5z,o|^ to the 
unperturbed field intensity, integrated over the extent of the tip above the slab. Field intensities in (b) and (e) are normalized 
to the maximum total field intensity |Eop. 



traces of the total -E-field energy in the cavity. 

Figure njc) displays a contour plot of the relative fre- 
quency shift Alu/lo as a function of the lateral position 
of a silicon tip of 125 nm diameter, placed at a height of 
30 nm above the photonic crystal slab. As expected, 
the introduction of high index material in the cavity 
mode causes a red shift of the mode frequency, depend- 
ing on the tip position. The maximum relative frequency 
shift of 0.13% amounts to about 2 nm (or 260 GHz) at 
A = 1500 nm and occurs when the tip is placed above 
the central defect. At a lateral distance of 200 nm from 
the central hole, the shift is already less than 10% of 
this maximum value, but as the tip approaches the two 
shifted holes, again a large tuning of up to 200 GHz can 
be achieved. Figure |21 displays the dependence of the 
cavity resonance on the separation between the tip apex 
and the photonic crystal slab (ztip), for the case where 
the tip is aligned with the cavity center. The filled circles 
show that as the tip is approached from afar {z^ip > 0), 
the frequency shift grows exponentially with a 1/e length 
of d = 50 nm. When the tip is put through the slab 
(•Ztip < 0), the increasing amount of dielectric material in- 
serted in the mode profile continuously detunes the reso- 
nance to the red, saturating at a shift of Aiojuj = —3.7%. 
The black curve in Fig. El displays the variation of the 
unperturbed cavity mode intensity |Eo|^. The excellent 
agreement between this curve and the filled symbols for 
(ztip > 0) indicates that the frequency shift maps |Eop. 
This direct correspondence of the frequency shift to the 
mode profile is also evident from the strong similarity of 
the lateral images in Figs.^^b) and (c). 

An exciting aspect of cavity tuning with a tip is that 
a large frequency shift can be achieved without inducing 
a considerable cavity loss. Open symbols in Fig. |21 show 



that within the accuracy of the simulations the cavity 
Q of 13000 has not been affected at separations down 
to about 60 nm while the cavity resonance has been 
tuned by 100 GHz. Even at a height of 30 nm where 
Aw=260 GHz, the cavity Q remains as high as 5000. 
If the tip is pushed closer to or into the central defect 
hole, however, the continuous increase of the frequency 
shift is accompanied by a strong reduction in Q. In this 
regime, the tip perturbs the mode profile and strong out- 
of-plane scattering occurs. When the tip extends com- 
pletely through the defect slab and the distribution of 
dielectric material becomes more symmetric, the Q re- 
covers to about 1000. 

To gain quantitative insight into the modification of 
the cavity resonance, we have taken advantage of the es- 
tablished theory of dielectric perturbations in microwave 
cavities 0| . Suppose the tip has a dielectric constant 
of Ep and is contained within a volume Vp. Maxwell's 
relations imply the exact expression, 

Auj i _- Iv, h - R-eiES • Ep]dr + ^ Re[Sp] • nda 
~ 2AQ ^ /^Re[E* -Dp + H* -BJdr 

(1) 

for the relative frequency shift (Aw/oj) and induced scat- 
tering loss {AQ~^). Here e(r) is the dielectric constant 
of the unperturbed system and Eg^p), and Ho^p) are the 
complex electric and magnetic fields of the unperturbed 
(perturbed) system. The integral in the denominator 
runs over some volume V enclosing the entire system. 
The induced scattering loss is proportional to the inte- 
grated flux of Sp = E* X Ho - E|5 X Hp through the outer 
surface 5V of V . We note that this dependence is very 
different to that of the frequency shift that is proportional 
to an integral over the perturbing volume Vp. 
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FIG. 2: Relative frequency shift (solid symbols, left hand axis) 
and inverse Q (open symbols, right hand axis) as a function of 
the separation Ztip between the surface of the photonic crystal 
and a silicon tip of 125 nm diameter aligned with the central 
defect. The grey area indicates the extent of the membrane. 
For Ztip < the tip extends into or through the defect. The 
dashed line indicates Ztip = 30 nm corresponding to the data 
in Figs, n^b-d). The black curve (scaled) shows the unper- 
turbed mode profile |Eo|^. The horizontal line indicates 1/Q 
in absence of the tip. 



FIG. 3: Relative frequency shift (solid symbols, left hand axis) 
and inverse Q (open symbols, right hand axis) as a function of 
the effective polarizability per length Oca/d — 3(ep — l)/(ep + 
2)7rrp, for cylindrical tips positioned 30 nm above the central 
defect in Fig. 1. Symbol shapes indicate various tip materials 
(Si, Ti02, Zr02, polystyrene, and Si02 in order of decreasing 
Ep as indicated). The frequency shift is proportional to Qch 
(solid curve). The loss 1/Q can be fitted (dashed curve) with 
the sum of the loss without the tip (horizontal dashed line) 
and a contribution proportional to a^g. 



A simple picture for the frequency shift emerges in the 
weakly perturbative regime where ^ < 1 and the Q 
is not degraded, in our case corresponding to the range 
Zfip > 30 nm. In this regime, the denominator in Eq. 
is independent of the perturbation and the frequency 
shift is set by the overlap integral of the perturbed and 
unperturbed fields within the tip volume Vp only. The 
field Ep in the tip is given by Ep = 3Eo/(ep + 2), where 
the proportionality constant takes into account the local 
field effects Next, we separate the exponential z- 

dependence of Eq = Eo(r||)e~^/^'' from its dependence 
on the lateral tip coordinate r||, and assume that the un- 
perturbed field Eo(r||) is constant over the small tip cross 
section. We then find 

Lu 2Kavmax[e(r)|Sop] ' ^ ' 

Here Vcav = / e(i")|-EopcZr/max[e(r)|i?ol^] is defined as 
the cavity mode volume 0, and ags is the effective 
polarizability of the tip 

cteff = 3^^^-— -^Kff, (3) 

equal to the electrostatic polarizability of a volume VcS 
of material with dielectric constant Cp [T^ . The expo- 
nential decay of the cavity mode (with decay length d) 
limits the volume of the tip that contributes to the po- 
larizability to Ves = irrpd for a tip of radius Vp [13. This 
analysis reproduces all the features of our FDTD results 
concerning the frequency shift. Firstly, it yields the ex- 
ponential decrease of Alu/lo as a function of the tip-slab 



separation znp (see Fig. Secondly, it confirms the 
direct correspondence between the frequency shift and 
the unperturbed mode profile (see Figs.^b,c)). Thirdly, 
Eq. lI2Jl states that the induced frequency shift is inversely 
proportional to the mode volume, a feature that we have 
confirmed by simulations of different cavity designs (e.g. 
that of Ref A]). Further confirmation is the quantita- 
tive agreement of the frequency shifts in Fig. [3 with the 
mode volume Vcav = 0.02 /im'^ determined from the un- 
perturbed mode profile of our cavity. Finally, and most 
significantly, Eq. ((2Jl predicts that the frequency shift is 
proportional to acff. To investigate this we performed 
FDTD simulations for tips of various radii and materi- 
als, fixed above the cavity center at 30 nm. The closed 
symbols and the solid line in Fig.|3demonstrate that the 
frequency shift is proportional to the effective polarizabil- 
ity of the tip. In other words, in the region ztip > 30 nm 
the tip acts as a polarizable object that is weakly cou- 
pled to the cavity field, much in the same manner as an 
atom that can be described by a classical dipole of the 
appropriate polarizability within linear dispersion the- 
ory 10, 11]. This is a remarkable finding because the tip 
is not much smaller than the central hole of the micro- 
cavity and its mode volume, so that the validity of dipole 
approximation is not a trivial matter. 

Next we turn our attention to the origin of losses in- 
duced by the tip. The power radiated by an object of 
polarizabilit y a pff can be modelled as a Rayleigh scatter- 
ing process [ig, resulting in a quadratic dependence on 
acff. Indeed the open symbols in Fig. |3| display exactly 
this dependence when the tip is at the center of the cav- 
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ity. The quadratic dependence of the loss compared to 
the linear dependence of the resonance frequency on aos 
makes it possible to tunc the frequency of high-Q mi- 
crocavities without incurring prohibitive losses. Intrigu- 
ingly, however, the induced loss is not directly related 
to the strength of the unperturbed cavity mode at the 
position of the tip. A comparison of Figs, ^b) and (d) 
shows that there are tip positions where the losses are 
high even though the electric field intensity is small (up- 
per right and lower left corners). It turns out, as shown 
in Fig. me), that the cavity field at these points has a 
large z-component {Ez.o). The magnitude of |-Ez,oP at 
these points amounts to about 5% of the maximum value 
of |EoP attained above the center of the cavity. Sub- 
stantial loss is induced by scattering of the field into the 
tip, which is known to be much more efficient for the 
z-component of the field than for the field components 
perpendicular to the tip 0|. FDTD simulations of the 
perturbed fields further confirm the polarization-specific 
enhanced coupling into the tip 20]. The disparity be- 
tween induced loss and the cavity mode profile under- 
lines a fundamental feature of Eq. contrary to the 
case of the frequency shift, the scattering loss does not 
depend solely on the fields within the perturbing volume 
only. Instead, it is the far-field interference of the fields 
encoded in the integrated flux Sp that is the decisive fac- 
tor. In this sense, there is a close relation between the 
perturbation of a photonic crystal microcavity by a tip 
and the tuning of holes around the cavity as used recently 
for the optimisation of the quality factor |^ ^ |E| - 

In conclusion, we have shown that sharp dielectric tips 
in the near-field can be used to tune photonic crystal 
microcavities over a large range without inducing pro- 
hibitive losses. Our analysis shows that the tuning range 
scales with the ratio of the tip polarizability to the mode 
volume. Hence, this scheme is expected to be effective 
for any electromagnetic resonance localized in a volume 
comparable to the perturbing polarizability, as is the case 
in microspheres, microdisks or micropillars [^IHI^' One 
advantage of this scheme is that the tuning process does 
not influence the optical properties of the emitters em- 
bedded in the cavity. Furthermore, the rapid progress in 
nanotechnology makes it feasible to integrate a tip with 
a microcavity in a device (iH. Compared to existing 
methods for tuning photonic crystals ^] , tip tuning can 
create and reverse the frequency shift more quickly and 
in contrast to ultrafast optical tuning [2^ , the frequency 
shift can be maintained indeflnitely. In other words, tip 
tuning does not depend on the material properties of the 
crystal. The time scale for tuning is only limited by me- 
chanical resonances of the tip which can reach the MHz 
regime jl^l • This opens the possibility of using a photonic 
crystal microcavity as an optical switch for integrated op- 
tics applications. 
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(Color) Panel (A): total E-field intensity of the cavity mode in a plane 30 nm above the slab (reproduced from manuscript, 

Fig. 1 (b)). Lines indicate the directions of the crosscuts along x in (C)-(E) and along y in (F), and circles mark the tip 
positions in (D)-(F) as labelled. Panel (B): cavity loss versus tip position (manuscript, Fig. 1(d)). Panels (C)-(F): Contour 
plots of the total E-field intensity |Eo|^/ max(|Eo|'^) (top half of each panel) and the contribution of the ^-component E-field 
intensity |-Ez,o|'^/ max(|Eoj^) (bottom half of each panel). Color scales are logarithmic (shown on the right of the panels). 
Field intensities are plotted for vertical plane cuts through the photonic crystal membrane along directions indicated by 
dotted lines in (A). Compared to the unperturbed profile (C), panel (D) shows that if the tip is above the central hole 
(position I in panel (A)) light does not couple into the tip, but scatters above and below the slab. Panel (E) shows that 
efficient coupling into the tip occurs when the tip is located in between two untuned holes (position II in panel (A)), leading 
to a large reduction in Q despite weak overall field intensity at position II. This coupling into the tip is selective for lobes of 
the z-component of the E-field, and so does not occur for position I amd III (panel F). In all cases. Si tips of 125 nm diameter 
(dotted outlines in (D)-(F)) were considered, at a height 30 nm above the crystal slab (top and bottom face indicated by 

dashes) . White bars in all panels indicate a = 420 nm. 



